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ABSTRACT 

We use high spatial resolution observations and numerical simulations to 
study the velocity distribution of solar photospheric magnetic bright points. The 
observations were obtained with the Rapid Oscillations in the Solar Atmosphere 
instrument at the Dunn Solar Telescope, while the numerical simulations were 
undertaken with the MURaM code for average magnetic fields of 200 G and 
400 G. We implemented an automated bright point detection and tracking al- 
gorithm on the dataset, and studied the subsequent velocity characteristics of 
over 6000 structures, finding an average velocity of approximately 1 kms~^, with 
maximum values of 7 kms~^. Furthermore, merging magnetic bright points were 
found to have considerably higher velocities, and significantly longer lifetimes, 
than isolated structures. By implementing a new and novel technique, we were 
able to estimate the background magnetic flux of our observational data, which 
is consistent with a field strength of 400 G. 



Subject headings: Sun: activity — Sun: evolution — Sun: photosphere 
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Introduction 



Granulation forms the instantly recognisable "patchwork" pattern of the solar photo- 
sphere and dominates quiet Sun regions. The dark intergranular lanes are formed as a result 
of convective downflows and are the regions where magnetic bright points (MBPs) can be 
observed. It is the horizontal transfer of magnetic flux from the centre of granules, and 
into the lanes, that forms MBPs. In the MBPs, the magnetic flux clumps tog ether to form 



small rnagnet ic concentrations with field strengths of the order of a kiloGauss (jStenfldll985 



Solankil 119931 ). These are some of the smallest features currently observable on the solar 
surface, and due to their dynamic natur e, can provide a con duit for the transfer of kinetic 
energy into the upper solar atmosphere (Ide Wijn et al.ll2009l ). 



Due to the increased computational power available in recent years, significant progress 
has been made in understanding the complex interplay between convective energy transport 
and radiation in these concentrated magnetic eleme nts. MBPs are best observed as strong 
intensity enhancements in G-band intensity images (IMuUer fc Roudierlll984j : iBerger fc Title 
19961) which have intensity peaks between 0.8 — 1.8 times the mean photospheric intensity 
( jSanchez Almeida et al. 2004 : Langhans et al. 2002). These enhancements are due to an 
increase in continuum intensity, caused by the continuum formation layer being depressed 
into the deeper, hotter layers of the solar photosph ere, where regions o f strong magnetic fie ld 
concentrations and reduced CH abundances exist (ISteiner et al.ll200ll : IShelyag et al.l 120041 ). 



The significance of MBPs on solar atmospheric energy tra nsport and dynamics ha s led to 
several studies on their area coverage and size distributi on. ISanchez Almeida et al.l (120041 ) 



found that MBPs have a minor axis of 135 km, while Wiehr et al 



predominant diameter of 160 ± 20 km. In a recent study, ICrockett et al. 



(2004) determined a 



(120101 ) utilised high 



spatial resolution observations and magneto-hydrodynamic (MHD) simulations to conclude 
that the MBP area distribution peaks at ^45000 km^, with a sharp decrease in occurrence 
for smaller surface areas. The area of the smallest MBPs is defined by the width of the inter- 
granular lanes, which is in turn limited by the balance of radiative and convective energies 
in the magneto-convective processes. 

MBP velo cities are induced by t he expansion of granules, and are typi cally in the range 



of 1 — 3 km s ^ ( IBerger fc Titlelll996l : IUtz et al.ll2010l ). It has been shown by lChoudhuri et al. 



( 119931 ) that magnetic footpoints with velocities greater than ~2 km s ^ can excite magneto- 
sonic kink waves, which can transport sufficient energy to heat the loca lised quiet corona . 



under a two layer atmospheric approximation. Al so, recent observ ations jBonet et al" 



2008 



Wedemeyer-Bohm fc Rouppe van der Voortl[2009[ ) and simulations (jShelyag et al.ll2011al ) sug 
gest the presence of vortex motions in the photosphere, with a poss ible connection to the 
spiral tracks of MBPs and chromospheric swirls. Very recent studies by lShelyag et al.l (l2011bl ) 
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and iFedun et al.l (1201 ll ) have indicated how vortex motions generate a significant amount of 
Poynting flux directed outwards from the photosphere, and, as a result, may be the source 
of various observed MHD wave modes. 

In this paper, we use high spatial and temporal resolution observations, in addition to 
numerical simulations, to determine the velocity distribution of a large sample of MBPs. The 
observations and numerical simulations are described in §|2l while the methodology used, and 
the values obtained for the velocities of MBP structures, are detailed in § [3J As our tracking 
algorithm can detect and monitor bright point chains, as well as isolated brightenings and 
merger events, we believe that this is a unique study of the dynamics of MBPs in the solar 
photosphere. Differences between the velocity characteristics of non-merging MBPs, and 
those that undergo mergers with other bright points, are discussed in § Ei Finally, our 
concluding remarks are given in § |H 



Observations and Numerical Simulations 



The data emplo yed in this study were obtained using the Rapid Oscillations in the Solar 
Atmosphere (ROSA: lJess et al.ll2010l ) instrument, which is installed as a common-user facility 
at the 76 cm Dunn Solar Telescope (DST), in New Mexico, USA. Observations were obtained 
during a period of excellent seeing on 2009 May 28, using a 9.2 A wide filter centred at 4305 A 
(G-band). We observed a 70" x 70" quiet Sun region at disk centre for ~50 minutes, achieving 
diffraction-l imited imaging with 0".069 pixel"^. The images were reconstructed using Speckle 
algorithms (jWoger et al.l 120081 ). while image de-stretching w as perforrned us i ng a 40 x 40 
grid (equating to a ~ 1".7 separation between spatial samples; I Jess et al.ll2007ll2008l ). These 
processes were implemented to remove the effects of atmospheric seeing from the dataset. 
G-band images were taken at a raw cadence of 0.033 s, while after speckle reconstruction the 
cadence was reduced to 0.528 s. Reconstructed images were then binned into consecutive 
groups of four to improve the signal-to- noise and reduce the overall volume of the dataset, 
providing a final image cadence of 2.1 s. 

Simulated G-b and images we r e pro duced using the detailed radiative transport tech- 
nique described by IShelyag et al.l ( 120041 ). with the solar photospheric magneto-convection 
model s for the radiative transport calculations provided by the MURaM radiative MHD 
code ( IVogler et al.l 120051 ). A computational domain of size 12 x 12 x 1.4 Mm^, was em- 
ployed for the simulations, resolved by 480 x 480 x 100 grid cells, providing a horizontal 
two-pixel resolution of 50 km. The level corresponding to the visible solar surface is located 
approximately 600 km below the upper boundary of the domain. Side boundaries of the 
domain are periodic, while the upper boundary is closed for vertical and stress-free hori- 
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zontal plasma motions and the bottom boundary is transparent. The nume rical solution is 
stabi l ised against numerical instabilities usin g hyperdiffusive source terms (ICaunt &: Korpi 
2OO1I : I Voder et aliboosl : Ishelvag et allboosi ). 



In the initial stage of the numerical simulation, we use a well-developed non-magnetic 
photospheric convection snapshot, where we introduce uniform unipolar vertical magnetic 
fields of 200 G and 400 G, representing two different levels of solar magnetic activity. Fol- 
lowing the injection of magnetic fields into the simulation, they quickly become horizontally 
transfered into the intergranular lanes, creating magnetic field concentrations with a strength 
of up to 2 kG at the photospheric level. After a few convective turnover timescales, we record 
219 and 177 magneto-convection snapshots, at 200 G and 400 G, respectively, for further 
analysis. To obtain the G-band intensities emerging from the simulated photospheric snap- 
shots, a spectral band in the wavelength range from 4295 A to 4315 A, which contains 241 
CH molecular lines and 87 atomic lines, was calculated for each vertical ray in the snapshots, 
and convolved with a filter function centred at 4305 A with a width of 9.2 A. 



3. Analysis and Results 

The MBP s anal ysed in this paper were detected using the algorithm developed by 



Crockett et al.l (120101 ). This automated algorithm uses the intergranular lanes to produce a 
binary map of bright features. A compass search is employed to disentangle features within 
the lanes, i.e. MBPs, from bright granules. Finally, intensity thresholding across the feature 
is used to determine more accurately the size and shape of the MBP. Long-lived features 
are then established, and stabilised, by examining subsequent frame s to observe the MBP's 



evolution in time. More details of these procedures may be found in lCrockett et al.l (120101 ). 



Figure [T] shows the full ROSA field-of-view, while Figure [2](a) shows an expanded version 
of the 16" X 16" region highlighted in Figure [H together with sample images of the 200 G 
(Fig. 12)3) and 400 G (Fig. [2t) simulations. The algorithm produces an outline of the centre 
of gravity of each MBP observed, as well as a direct measurement of the area of each feature. 
A 30 s cut-off was employed to remove those MBPs which are short-lived and unlikely to 
have a significant horizontal velocity component, as well as assisting the removal of any false 
detections that may adversely affect the results. Figure ?? shows snapshots from two videos 
that have been included as an online resource. The first depicts the tracking capabilities of 
the algorithm for an individual MBP throughout its 336 s lifetime, while the second displays 
how an MBP from the 200 G MHD simulation is tracked. 



Our observed G-band dataset included a total of 6236 MBPs with hfetimes exceeding 
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30 s. These MBPs were found to have an average velocity of approximately 1 kms~^, and 
a mean lifetime of 91 s. The highest velocity recorded in this study was 7 kms~^, a value 
comparable to the photospheric sound speed, while the longest MBP lifetime observed was 
approximately 20 minutes. Similar values were obtained when the same methodology was 
applied to the numerical simulations. A total of 721 MBPs were detected in the 16" x 16" 
200 G simulation, with an average velocity of 0.6 kms^^, a maximum of 4.1 kms~^, an average 
lifetime of 108.2 s, and a maximum of just over 11 minutes. In the 400 G simulation, we 
detected 1024 MBPs, with an average velocity of 0.65 kms~^ and a maximum of 6.2 kms^^, 
an average lifetime of 168.4 s and a maximum of nearly 17 minutes. Histograms indicating the 
distributions of velocities for all three data sets are shown in Figure [2t^d-f), with the results 
summarised in Table 1. A dependence of the velocity distribution on the average magnetic 
field strength (and thus on the solar magnetic activity level) suggests an alternative, and 
novel, way of estimating the net magnetic flux in the solar photosphere. Hence, as the best 
agreement with our observations is found for the 400 G simulation, we can infer the average 
magnetic flux of our observations as being close to 400 G. Small discrepancies between 
our observed and simulated datasets may be due to the large MBP sample present in our 
observations (6236 versus 1024). 

The distribution of velocities, for both the observed and simulated datasets, are shown 
as histograms in Figure [2t^d-f). These demonstrate that the majority of MBPs have velocities 
between and 1 kms~^, while some MBPs (i.e. 30% in the observations, 6% in the 200 G 
simulated data and 6% in the 400 G) have velocities in excess of 2 kms~^. The considerable 
fraction of MBPs with velocities greater than 2 kms~^ is important, as rapid bursts in their 



motion may induce magnetosonic kink waves (jChoudhuri et al.lll993l ). As shown previously. 



these kink waves may act as a conduit for imparting energy into the upper solar atmosphere. 

Although there is considerable agreement between the results obtained in the observed 
and simulated data, there are some slight discrepancies that need to be addressed for a 
more complete picture. For example, the maximum velocity, the maximum lifetime and the 
proportion of MBPs with velocities greater than 2 kms~^ are all smaller in the simulated 
datasets. A possible reason is that the net (or average) magnetic field applied within the 
computational domain remains constant throughout the 200 G or 400 G time series, whereas 
the flux present in the observed data may constantly change with time. This dynamic nature 
of flux generation and annihilation could directly lead to the higher velocity distributions 
detected for the observed MBPs. 



An alternative algorithm based on local correlation tracking (LCT; IWelsch et al.l 12004 



Matloch et al.ll2010l ) was applied to the observations and simulated time series, and yielded 
similar velocity characteristics to within 0.3 kms~^. A significant number of MBPs, ~21% 
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of those detected in the observations, were involved with a merging event. Only a few of 
mergers (31) are observed in the 200 G simulation, whereas 207 are detected in the 400 G 
simulation. Comparisons with an active region area would be a useful study in the future 
to observe how the frequency of merging events in such a feature compares to the quiet Sun 
region discussed here. 

There is only a slight difference between the velocities obtained for merging events and 
those of isolated MBPs. The average velocities evaluated for all the merging events in the 
observations was found to be 1.6 kms"^, slightly above the measured value for isolated 
MBPs. However, the average lifetime of the merging MBPs is 201.2 s, more than twice the 
value for non-merging MBPs. These results are echoed in the values for the 200 G simulated 
dataset, although only 31 mergers were witnessed out of the 721 detections. Here, the average 
velocity of the merging MBPs is 0.9 kms~^, with an average lifetime of 227 s. For the 400 G 
simulation time series, the average velocity was found to be 0.8 kms~^, with an average 
lifetime of a merging MBP of 329 s (determined over a total of 207 merging events). These 
results show that there is little difference in the velocity of a merging MBP when compared 
with its non- merging counterpart. However, the former do appear to have longer lifetimes. 
This effect may arise due to an increase in stability of the MBPs when they undergo a merger. 
One possibility to account for the marginal elevation of velocity characteristics for merging 
MBPs is the barycenter method used to estabhsh the merging velocity. During a merger, the 
shift of the centre of gravity as the MBPs unite may lead to a larger recorded velocity value. 
This may also explain why we see such a high maximum velocity in the observed data set. 
It is most hkely an artefact of MBP geometry changes, which induce a considerable motion 
on the barycenter of the MBP structure. Thus, some velocities may be overestimated by 
this method. However, employing LCT methods on the same MBP groups indicate that a 
fraction of these structures do in fact demonstrate velocities comparable to the sound speed. 



4. Concluding Remarks 

We have used high resolution observations, and state-of-the-art MHD simulations, to 
determine the velocity distribution of magnetic bright points (MBPs) on the solar surface. 
The average velocity established from our observations is 1 kms~^, averaged over a total of 
6236 MBP structures, with a typical lifetime of 91 s. A total of 721 MBPs were tracked in 
the 200 G numerical simulation, with an average velocity of 0.6 kms^^ and a mean lifetime 
of 108.2 s. The slight differences between the observed and 200 G simulated datasets can be 
attributed to the initial value of magnetic field applied to the MHD simulation. Using an 
initial magnetic field of 400 G, the subsequent numerical simulation yielded values closer to 
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the observational distribution. The average velocity and average lifetime of the 1024 MBPs 

detected in the 400 G data is 0.65 kms~^ and 168.4 s, respectively. When local correlation 
tracking (LCT) was applied to each data set, the difference in MBP velocities was found to 
be ~0.3 kms~^, showing a level of consistency between the two independent techniques. 

For the merging events witnessed in the observational dataset, the average velocity was 
found to be 1.6 kms"-*^, with a mean lifetime of 201.2 s. The 200 G simulated dataset only 
produced 31 merger events, with an average velocity of 0.9 kms~^ and mean lifetime of 227 s. 
By contrast, the 400 G simulation had 207 mergers, with an average velocity of 0.8 kms~^ 
and mean lifetime of 329 s. All of the data sets show that merging events exhibit longer 
average lifetimes, which we believe is a result of greater MBP stability as a result of the 
merger. Although the 200 G simulation produce similar MBP lifetimes, our data are most 
consistent with the 400 G simulation in both appearance and in MBP merger lifetimes. We 
intend to expand on this study by observing the velocity characteristics of MBPs and merger 
events in the vicinity of an active region. 

Importantly, we have presented a novel technique for estimating the average photo- 
spheric magnetic flux through comparison of observations and MHD simulations. By scaling 
the initial magnetic flux injected into the numerical simulation to replicate the MBP dy- 
namics established in our observational time scries, we can directly estimate the average 
background magnetic flux in the solar feature studied, without the necessity of vector mag- 
nctograms. Following this methodology, we find our observational data is consistent with an 
average magnetic flux of ~400 G. 
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the award of a Post-Doctoral Fellowship. Effort sponsored by the Air Force Office of Scientific 
Research, Air Force Material Command, USAF under grant number FA8655-09-13085 

REFERENCES 

Berger, T. E., & Title, A. M., 1996, ApJ, 463, 365 

Bonet, J. A., Marquez, I., Sanchez Almeida, J., Cabello, I., & Domingo, V. 2008, ApJ, 687, 
L131 



-8- 

Caunt, S. E., & Korpi, M. J., 2001, A&A, 369, 706 

Choudhuri, A. R., Dikpati, M., & Banerjee, D., 1993, ApJ, 413, 811 

Crockett, P. J., Mathioudakis, M., Jess, D. B., Shelyag, S., Keenan, F. P., & Christian, D. J., 
2010, ApJL, 722, L188 

de Wijn, A. C, Stenflo, J. O., Solanki, S. K., & Tsuenta, S., 2009, Space Sci Rev., 144, 275 

Fedun, V., Shelyag, S., Verth, G., Mathioudakis, M., & Erdelyi, R. 2011, Annales Geophys- 
icae, 29, 1029 

Jess, D. B., Andic, A., Mathioudakis, M., Bloomfield, D. S., & Keenan, F. P. 2007, A&A, 
473, 943 

Jess, D. B., Mathioudakis, M., Christian, D. J., Keenan, F. P., Ryans, R. S. I., & Crockett, 
P. J., 2010, Sol. Phys, 261, 363 

Jess, D. B., Mathioudakis, M., Crockett, P. J., & Keenan, F. P., 2008, ApJL, 688, L119 

Langhans, K., Schmidt, W., & Tritschler, A., 2002, A&A, 394, 1069 

Matloch, L., Cameron, R., Shelyag, S., Schmitt, D., & Schiissler, M., 2010, A&A, 519, A52 
MuUer, R., & Roudier, T., 1984, Sol. Phys., 94, 33 

Sanchez Almeida, J., Marquez, I., Bonet, J. A., Dommguez Cerdena, I., & MuUer, R., 2004, 
ApJL, 609, L91 

Shelyag, S., Schiissler, M., Solanki, S. K., Berdyugina, S. V., Vogler, A., 2004, A&A, 427, 
335 

Shelyag, S., Fedun, V., & Erdelyi, R. 2008, A&A, 486, 655 

Shelyag, S., Keys, P. H., Mathioudakis, M., & Keenan, F. P., 2011a, A&A, 526, A5 

Shelyag, S., Fedun, V., Keenan, F. P., Erdelyi, R., & Mathioudakis, M. 2011b, Annales 
Geophysicae, 29, 883 

Solanki, S. K., 1993, Space Sci. Rev., 63, 1 

Steiner, O., Hauschildt, P. H., & Bruls, J., 2001, A&A, 372,L13 

Stenflo, J. O., 1985, Sol. Phys., 100, 189 



-9- 

Utz, D., Hanslmeier, A., MuUer, R., Veronig, A., Rybak, J., & Muthsam, H., 2010, A&A, 
511, A39 

Vogler, A., Shelyag, S., Schiissler, M., Cattaneo, F., Emonet, T., & Linde, T., 2005, A&A, 
429, 335 

Wedemeyer-Bohm, S., & Rouppe van der Voort, L. 2009, A&A, 507, L9 
Welsch, B. T., Fisher, G. H., Abbett, W. R, & Regnier, S. 2004, ApJ, 610, 1148 
Wiehr, E., Bovelet, B., & Hirzberger, J., 2004, A&A, 422, L63 
Woger, F., von der Liihe, O., & Reardon, K., 2008, A&A, 488, 375 



This preprint was prepared with the A AS macros v5.2. 



-20 



-10 10 

X (arcsecs) 



20 



Fig. 1. — The full 70" x 70" ROSA field-of-view of the observational G-band data used in this 
study. The dashed white box indicates a 16" x 16" region which is expanded in Figure EJ^a). 
The axes are in heliocentric arcseconds, with (0,0) representing the centre of the solar disk. 
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Fig. 2. — The upper panels contain 16" x 16" snapshots of the three datasets used in this 
study, (a) is an expanded view of the boxed region in Figure [T] while (b) and (c) represent 
the 200 G and 400 G simulations, respectively. The velocity distributions for each dataset 
are shown below their corresponding snapshot in panels (d-f). 
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Tal)lo 1: Suiiiiiiary of MBP cliaractoristics. 



Data Set 


ROSA Observations 


200G Sims 


400G Sims 


FOV size 


70" X 70" 


16" X 16" 


16" X 16" 


Average cadence (s) 


2.112 


8.7 


6.7 


Number of MBPs detected 


6236 


721 


1024 


Mean lifetime (s) 


91 


108 


168 


Average velocity (kms^^) 


1.0 


0.6 


0.65 


Maximum velocity (kms^^) 


7.0 


4.1 


6.2 


Percentage of detections that merge 


21% 


4% 


20% 


Average velocity of mergers (kms"-*^) 


1.6 


0.9 


0.8 


Average lifetime of mergers (s) 


201 


227 


329 



